Constraining Anisotropic Baryon Oscillations 
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We present an analysis of anisotropic baryon acoustic oscillations and elucidate how a mis-estimation of the 
cosmology, which leads to incorrect values of the angular diameter distance, Úa, and Hubble parameter, H, 
manifest themselves in changes to the monopole and quadrupole power spectrum of biased tracers of the den- 
sity field. Previous work has focused on the monopole power spectrum, and shown that the isotropic "dilation" 
combination S A H~ X is robustly constrained by an overall shift in the scale of the baryon feature. We extend 
this by demonstrating that the quadrupole power spectrum is sensitive to an anisotropic "warping" mode aU-ff, 
allowing one to break the degeneracy between úa and H . We describe a method for measuring this warping, ex- 
plicitly marginalizing over the form of redshift space distortions. We verify this method on N-body simulations 
and estimate that ^aH can be measured with a fractional accuracy of ~ (3/W)% wherethe survey volume is 
estimated in /i _3 Gpc 3 . 

PACS numbers: 



I. INTRODUCTION 

Oscillations in the baryon-photon fluid prior to decoupling 
leave an imprint in both the cosmic microwave background 
angular power spectrum and the power spectrum of the mat- 
ter with a characteristic length scale fljj] . Measurements of 
this scale at low redshifts, for example in large galaxy red- 
shift surveys, provide a standard ruler for constraining the 
expansion history of the Universe. A unique feature of this 
method is that it can constrain both the angular diameter dis- 
tance Úa(z) and the Hubble parameter H(z) at a given red- 
shift, z. However, the majorityof current work, both theo- 
retical and observational Ü H |1 H H, H, has focused 
on angle-averaged baryon oscillations which only constrain 
the combination d\H~ l . While there have been initial ef- 
forts to use the full two-dimensional information [fl~oh - these 
have assumed knowledge of the form of redshift-space dis- 
tortions lílll Hlh . An open question is how to use the full 
two-dimensional information robustly, i.e. without assuming 
knowledge about the particular form of redshift space distor- 
tions. 

This note attempts to fill this gap. We recast the rescaling 
of the transverse and line-of-sight directions as "dilation" and 
"warping" modes, and demonstrate that a Legendre expansion 
of the two-dimensional power spectrum (or correlation func- 
tion) can effectively isolate these two modes. In particular, 
we show that the warping mode can be constrained by match- 
ing to a template determined by the angle-averaged power 
spectrum, without any assumptions about redshift space dis- 
tortions. 



H. DILATIONS VS WARPING 

Theory: In order to understand the effects of incorrect angu- 
lar and radial distance scales on the baryon oscillation signal, 
we start with the two dimensional power spectrum P(fc,/z) 
where fc is the amplitude of the wave-vector with components 
k± and k\ | perpendicular and parallel to the line of sight, and 



/i = feii/fcis the cosine of the angle between the line of sight 
direction and k. We use the flat-sky approximation throughout 
for simplicity. and work in Fourier space, although analogous 
expressions exist configuration space. It is also convenient to 
expand the angle dependence of P(k, ¡jl) in Legendre polyno- 
mials Ve(fi), 
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while e constrains the distance ratio (k±/k\ 
(d^Tí) -1 / 3 . This notation separates isotropic a deformations 
("dilation") from anisotropic e deformations ("warping"). The 
measurable effect of a dilation is to scale fc by a, or slide the 
Pg along the fc— axis. Since baryon oscillations imprint a fea- 
ture with a known scale on the í = multipole, it is straight- 
forward to determine a once this feature has been detected. 
Indeed, this is how baryon oscillations have been used most 
commonly to date as a standard ruler. Furthermore, since one 
is matching a feature in the power spectrum (or correlation 
function), this probe is robust to a large number of potential 
systematics. However, it should be emphasized that a does 
not constrain d^ and H independently, but rather the combi- 
nation d A H~ x . A dilation also changes the normalization of 
the power spectrum, but this is perfectly degenerate with the 
bias of the mass tracer. 

In order to understand the effect of e, we consider a toy 
model where P¿>2 = 0, a reasonable approximation on large 
scales. For small deformations e< lwe obtain, to first order 
in e, 
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Substituting into Eq. — we obtain 
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where we ignore any hexadecapole terms. In the region of 
the baryon oscillations, the two corrections to the £ = 
power have about equal magnitudes but opposite sign, result- 
ing in an almost exact cancellation, and therefore suppress- 
ing the correction term. As a result, the monopole power is 
effectively unchanged, as one might have naively expected; 
Fig. Q] demonstrates this explicitly for a suite of simulations 
described below. 

There are three effects of warping on the £ = 2 component. 
The first is a scaling of its amplitude. However, as a practical 
matter, the change in the amplitude will be degenerate with 
uncertainties in our understanding of redshift-space distor- 
tions and therefore must be marginalized over. The g?P 2 /d ln k 
term can be recast as a scale shift in the power spectrum, sim- 
ilar to the more familiar scale shift from the dilation described 
above. In principle, the existence of a feature in the power 
spectrum allows one to constrain such a shift, even marginal- 
izing over uncertainties in the broad shape of the galaxy power 
spectrum. However, Eq. — shows that the size of the shift is 
suppressed by a factor of ~ 2. This, combined with the fact 
that the amplitude of the feature is suppressed in the £ = 2 
component and the larger errors for 1=2 reduces the strength 
of statistical power of this term (we revisit this more quan- 
titatively below). The final term mixes in the derivative of 
the £ = component. Normally, such a mixing would have 
smoothly distorted the shape of the £ = 2 component, which 
the uncertainties in redshift space distortions would have then 
forced us to marginalize over. However, the oscillations im- 
print a distinctive feature in the derivative term that one can 
robustly search for. Note that the above three effects are noth- 
ing but parts of the Alcock-Paczynski lf~1 effect; the only 
difference with its more classical application being that the 
presence of a feature breaks the traditional degeneracies with 
broad-band shape and amplitude. 

In order to demonstrate that our toy model captures the es- 
sential phenomenology, we consider a suite of 40 particle- 
mesh N-body simulations each of 1024 3 particles with a 
2048 3 force mesh, run assuming a ACDM cosmology with 
Ü M = 0.25, fl A = 0.75, Q b h 2 = 0.0224, h = 0.7, n = 0.97 
and crg = 0.8 in a periodic 2 /i _1 Gpc box; the precise details 
of these simulations will be presented elsewhere |15]. For 
definiteness we focus on the output at z = 0.3, although this 
choice was arbitrary. We begin by considering the dark matter 
clustering. For each simulation, we warp the dimensions of 
the volume as described above, before computing the £ = 0,2 
power spectra (and correlation functions). The averages over 
the 40 boxes are shown in Figs. Q] and — In order to esti- 
mate the terms in Eq. — we Savitzky-Golay filter fj~ the un- 
warped power spectra with a fourth order, nine-point filter to 
estimate both the smoothed power spectra and its derivatives. 
The predictions for l = 2 are also plotted in Fig.[_ The model 
clearly captures the behavior of the £ = 2 power spectrum 
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FIG. 1: The effect of warping on the 1 = and 2 Legendre mo- 
ments of the power spectrum, measured by averaging over the 40 
PM simulations discussed in the text. To better highlight the baryon 
oscillations, we have divided out by the "no-wiggle" power spectrum 
of li~ . The filled squares [black] are the unwarped power spectrum, 
the triangles [blue] are warped by 2%, the squares [red] by 5%, and 
circles [green] by 10%. The solid lines are the predictions for £ = 2 
of the simple model in Eq. — The predictions for £ = are equally 
accurate, although we don't plot them to avoid cluttering the plot. 
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FIG. 2: As in Fig.fT) but for the redshift-space correlation function. 
Note that we assume £(r, /i) = £oVo — £,%Pi + so that the multi- 
poles are mostly positive. Note that the feature in £2 appears at larger 
separations than the comparable feature in £o- 



3 





Dark Matter 


Galaxies 


Warp 


n 


^max 


<c> 


Xdof 


(e) a £ 


2 

Xdof 





3 


0.2 


0.0 1.0 


0.93 


0.1 0.9 


0.91 


1 


3 


0.2 


0.9 1.0 


0.94 


1.1 1.0 


0.92 


2 


3 


0.2 


2.1 1.0 


0.97 


2.3 1.0 


1.01 


3 


3 


0.2 


3.2 1.1 


0.94 


3.2 1.0 


1.00 


5 


3 


0.2 


5.4 1.3 


0.96 


5.2 1.0 


1.08 


10 


3 


0.2 


9.9 1.6 


1.29 


9.3 1.1 


1.82 


1 


3 


0.1 


0.9 2.1 


1.03 


1.2 2.2 


1.13 


1 


3 


0.2 


0.9 1.0 


0.94 


1.1 1.0 


0.92 


1 


3 


0.3 


0.9 1.0 


0.69 


1.8 0.7 


0.69 


1 





0.2 


0.9 1.0 


0.89 


-3.9 0.6 


1.99 


1 


1 


0.2 


1.0 1.0 


0.92 


0.5 2.0 


0.94 


1 


2 


0.2 


0.9 1.0 


0.94 


1.1 1.0 


0.92 


1 


3 


0.2 


0.9 1.1 


1.02 


1.1 1.0 


0.98 



TABLE I: Constraining the warping e of the simulation boxes using 
the method outlined in the text for the dark matter and galaxies. All 
the warps are reported in per cent. The first group of numbers demon- 
strate that we correctly recover the warping, except for large warps 
where the model is known to be no longer accurate; a simple iter- 
ative procedure would work in such cases. The remaining numbers 
demonstrate that these results are insensitive to the particular choices 
of the bias and k cut-off, except for extremely restrictive cases. 



when warped and is a good fit for small warping although, not 
surprisingly, it starts to deviate as the warping increases. 
Implementation: The above discussion suggests a procedure 
to constrain the warping mode. Since the £ = mode is 
mostly unaffected by warping, we assume that we use it to 
correct for the isotropic shifts. Motivated by Eqs.|U we then 
fit P2{k) with a function of the form, 
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where the hats denote undistorted quantities. We marginalize 
over smooth shape distortions with a general bias function, 



B n (k) 
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Since the £ = power spectrum is mostly unaffected by warp- 
ing, we determine Po directly from the actual warped mea- 
surements even though Eq.[5]refers to the unwarped spectrum; 
the differences between the two are suppressed at linear or- 
der in e. For l = 2, we suggest that P 2 be determined from 
simulations, and allow B n (k) to correct for any template mis- 
matches. In what follows, we determine P 2 and the derivative 
of Po by Savitzky-Golay filtering the unwarped average l = 2 
and warped l = power spectra respectively; the parameters 
of the filter are the same as used previously. 

Table|I]shows the results of fitting Eq.[5]to the dark matter 
power spectrum in the simulations described earlier. The al- 
gorithm correctly recovers the input warp, except as expected 
for the case of large warps, where Eqs.[4]and [5]are no longer 
valid. For such cases, one would simply iterate the above pro- 
cedure. We also note that for such large warps, the £ = 
spectrum is also significantly distorted and is therefore eas- 
ily diagnosed. Table [1] also demonstrates that the recovered 



parameters are insensitive to either the form of B n or the cut- 
off in k, except when we exclude significant amounts of the 
baryon oscillation signal (fc max = 0.1/iMpc -1 ). This sug- 
gests that our constraints are from the features in the power 
spectrum, and not the broad-band shape. 

We separate the constraints from the shift in P 2 from the 
mixing of dPo /dhik by fitting for each separately in Eq. [5] 
instead of fixing the relationship between them. The shift is 
poorly constrained with an error of ~ 8%, while the error on 
the amplitude of the mixing term degrades to~ 1.2-1.3% 
depending on the particular case we consider. The two are 
significantly anti-correlated, as might be expected from Eq.[4] 
explaining the non-trivial improvement when they are com- 
bined. 

In order to test our algorithm on a galaxy-like sample, we 
populate the simulations with mock galaxies using a local- 
density dependent prescription designed to match the number 
density and clustering of a luminous red galaxy sample fl7Íl . 
Generalizing the above algorithm for galaxies is trivial, but for 
one subtlety - the choice of the template power spectra. While 
Po can be determined by the £ = warped galaxy spectrum, 
the form of P 2 must be determined from theory. We start with 
the £ = 2 dark matter power spectrum as above, but suppress 
it with a Gaussian e~( fc -R) to account for the effect of virial 
motions within halos. We do not simultaneously fit for the 
warp and R, but fix R to 5/i -1 Mpc appropriate to our partic- 
ular sample of galaxies, and allow B n (k) to correct for errors 
in the suppression term. In practice, one would determine the 
form of this suppression by fitting to the small-scale galaxy 
correlation function. The results are summarized in Table [1] 
as with the dark matter, we see that our algorithm correctly 
recovers the input warps, except in two cases. The method 
fails for /c max = 0.3/iMpc -1 , since (for this particular galaxy 
sample) smaller scales are significantly affected by virial mo- 
tions and Eq. [5] is no longer a good description of the power 
spectrum. In principle, this could be alleviated by e.g. collaps- 
ing clusters to remove the virial motions, or down-weighting 
these scales. Similarly, the method fails if one chooses too 
restrictive a form for B n (k); this is seen most spectacularly 
for n = 0, which assumes no scale dependence to the bias, 
although the best fit is disfavored in this case at greater than 
99.9%. Increasing n quickly corrects this problem, and we 
correctly recover the input warps. Both these cases emphasize 
the need to choose a flexible form to describe the connection 
between galaxies and dark matter. We have demonstrated that 
a simple choice works in the case considered here, but more 
comprehensive investigations must be the subjects of future 
work. 

Comparison to Fisher matrix estimates: A Fisher matrix cal- 
culation lfl8ll predicts that measuring baryon oscillations in a 
8 /i _3 Gpc 3 volume yields fractional errors on dA and H of 
0.7% and 1.8% respectively with a cross-correlation coeffi- 
cient of 0.4. This implies errors on a and e of 0.6% and 0.7% 
with a cross-correlation coefficient of 0.4. The Fisher error on 
e is 40% better than our estimates. This is not too surpris- 
ing, since the Fisher matrix calculation of IU8I1 makes optimal 
use of the entire power spectrum assuming a perfect param- 
eterized model (and also requires a subtraction of informa- 
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tion from broad-band power). The degradation of dark energy 
constraints will be smaller than this, since the dominant con- 
straint comes from the dilation mode; the precise degradation 
will depend on the redshift range of the experiment and the 
particular priors used. 

III. SUMMARY 

The standard ruler method for constraining the expansion 
history of the Universe, as exemplified by baryon acoustic os- 
cillations, can measure both cLa{z) and H{z) with high ac- 
curacy. However, the most commonly used angle-averaged 
measures constrain only the combination d^H^ 1 . We have 
shown that by measuring both the monopole and quadrupole 
power spectra P0.2, it is possible to break this degeneracy and 
measure cLa and H independently without assuming a form 
for redshift space distortions. 

By modeling the manner in which "warping" modes mod- 
ify P2, we have given an explicit method for measuring 
ÚaH, allowing us to break the degeneracy between Úa and 
H from the angle-averaged measurement. We have demon- 
strated this method on N-body simulations and estimate that 
the combination ÚaH can be constrained with an accuracy 
of ~ {3/W)% where the survey volume is estimated in 
/i _3 Gpc . Since our method relies only on the information 



contained in the acoustic feature, it is insensitive to the de- 
tails of the complex relation between mass and galaxies or the 
specific form of the redshift space distortions. This degrades 
the constraint on d^H by ~ 40% compared with traditional 
Fisher matrix estimates which assume perfect knowledge of 
redshift space distortions (even when they marginalize over 
the shape of the power spectrum). This is, of course, a con- 
servative assessment, and improvements in our abilities to ro- 
bustly the relationship between galaxies and matter will im- 
prove constraints. 
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